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ABSTRACT 

We argue that gamma-ray somxes observed in the direction of the Cyg 0B2 association 
in the GeV and TeV energy range are due to a pulsar which was created by a supernova 
a few tens of thousands years ago. The GeV emission is produced by a middle-aged 
pulsar, a factor of two older than the Vela pulsar. The TeV emission is produced by 
high energy hadrons and/or leptons accelerated in the pulsar wind nebula. We suggest, 
moreover, that the excess of cosmic rays at ^ 10^^ eV, observed from the direction 
of the Cygnus region, can also be related to the appearance of this very energetic 
pulsar in the Cyg 0B2 association. A part of relativistic hadrons, captured in strong 
magnetic fields of a high density region of the Cyg 0B2, produce neutrons and 7-rays 
in collisions with the matter. These neutrons can arrive from the Cyg 0B2 creating 
an excess of cosmic rays. 

Key words: gamma-rays: individual: Cyg 0B2 - ISM: clouds - cosmic rays ~ pulsars: 
general - gamma-rays: theory - radiation mechanisms: non-thermal 



1 INTRODUCTION 

The Cygnus OB2 is the largest massive galactic association 
at a relatively small distance of 1.7 kpc. It has the diam- 
eter of ~ 60 pc and the core radius of ~ 14 ± 2 pc. The 
total H2 mass of the Cyg OB2 derived from the CO survey 
is ~ 3.3 X 10^ M0 (Butt et al. 2003), and the total stellar 
mass is (4 - 10) x 10* M© (Knodlseder 2000). It contains 
many massive stars of the OB type (~ 2600 ± 400) and O 
type (~ 120 ± 20) (Knodlseder 2000). The typical main se- 
quence evolution lifetime of the massive O type stars is of 
the order of a few 10^ yrs and B type stars a few tens of 
10^ yrs. Massive O type stars should pass through the Wolf- 
Rayet phase and explode as supernovae. It is expected that 
very fast pulsars are created in explosions of such massive 
stars at a rate of about one per a few 10"* yrs within the 
Cyg 0B2. This rate can be much higher if most of these 
massive stars are of similar age. The pulsars, created ~ 10* 
yrs ago, should resemble the Vela type pulsars (e.g. Vela 
pulsar, PSR 1706-44), which are strong sources of pulsed 
7-rays in the hundred MeV-GeV energies (EGRET obser- 
vations, Kanbach et al. 1994, Thompson et al. 1996). These 
pulsars are also surrounded by relatively small pulsar wind 
nebulae (PWNe), which are sources of TeV 7-rays (Chad- 
wick et al. 2000, Kifune et al. 1995, Chadwick et al. 1998). 

The Cygnus region is difficult to observe in high en- 
ergies due to its very complex structure and location on 
the galactic plane. The analysis of the EGRET data from 
different viewing periods revealed a few sources in a rel- 
atively small field of view containing the Cyg OB2: 2EG 
J2033-h4112 (Mori et al. 1997), GeV J2035+4214 (Lamb 



& Macomb 1997), GRO J2034-f4203 (Reimer, Dingus & 
Nolan 1997), 3EG J2033+4118 (Hartman et al. 1999). The 
first source is interpreted by the authors as caused by the 
massive binary Cyg X-3. The other sources may be due 
to the same object in the Cyg OB2, according to Butt et 
al. (2003). 

Recently, the HEGRA group has reported the discovery 
of an unidentified, steady, TeV source at the edge of the 95% 
error circle of the source 3EG J2033-I-4118 (Aharonian et 
al. 2002) . The spectrum of this source is rather flat (differen- 
tial spectral index — 1.9±0.3atat ±0.3sys), with the flux about 
two orders below the flux observed from 3EG J2033-(-4118. 
Moreover, there is an indication that the TeV source is ex- 
tended with a radius of 5.6' (Aharonian et al. 2002), con- 
firmed in the new analysis of the HEGRA data by Rowell & 
Horns (2002). 

The discovery of TeV source has led to a more care- 
ful analysis of the X-ray observations of Cyg OB 2 from 
the ROSAT and Chandra telescopes (Mukherjee et al. 2003, 
Butt et al. 2003). None of the several X-ray sources detected 
by the Chandra in the region with ~ 11' diameter around the 
TeV source was particularly prominent (Butt et al. 2003). 
Mukherjee et al. (2003) notes only one interesting variable 
source which lies 7' from the TeV centroid with the flux of 
2 X 10~^^ ergs cm~'^ s~^. It might be a quiescent low-mass 
X-ray binary containing a neutron star or a black hole. How- 
ever, weak diffuse X-ray emission has also been observed, 
with the intensity of 1.3 x 10~^^ ergs cm~^ s~^ in 0.5-2.5 
keV and 3.6 x 10"^^ ergs cm"^ s"^ in 2.5-10 keV (Butt et 
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al. 2003). This is more than ten times as bright as the sum 
of all the point-like sources. 

Recent analysis of the arrival directions of cosmic rays 
with energies of ^ 10^* eV by the AGASA group shows 
the excess of particles from the direction of the Galactic 
Genter (Hayashida et al. 1999a, confirmed by the SUGAR 
analysis, Bellido et al. 2001), and from the direction of the 
Cygnus region. The second excess has a significance of 3.9 
cr (3401 events with the background of 3148) in the energy 
range 10^* - 10^** "* eV (Hayashida et al. 1999b). We suggest 
that this excess of EHE particles can also be related to the 
high energy processes occurring in the massive Cyg 0B2 
association. 

The possibility that 7-rays are produced by massive 
young stellar clusters, as a result of acceleration of parti- 
cles in the winds or shocks produced by the early type stars, 
has been frequently discussed in the past (with the applica- 
tion to the Cyg OB2, see e.g. Chen, White & Bertsch 1992, 
Manchanda et al. 1996, Benaglia et al. 2001). Aharonian 
et al. (2002) discusses such general scenario, suggesting that 
the TeV emission may originate in collisions of hadrons with 
a dense cloud. These authors also propose that 7-ray emis- 
sion may originate in the inverse Gompton scattering by 
electrons, which are accelerated by a jet-driven termination 
shock expelled from a microquasar, e.g. Gyg X-3 (see Aharo- 
nian & Atoyan 1998). Butt et al. (2003) argue that the TeV 
emission is more likely to originate in collisions of hadrons 
with the matter than in the leptonic processes. Hadrons 
might be accelerated to PeV energies in collisions of strong 
stellar winds. 

In this paper we interpret the variety of high energy pro- 
cesses connected with the Cyg 0B2 association by applying 
the model in which the basic role is played by a young pul- 
sar, born in this association about ~ 10* yrs ago. It is shown 
that the observed GeV emission can be due to the Vela type 
pulsar present inside or close to the Cyg OB2 (Sect. 3). The 
TeV 7-ray emission originates in the pulsar wind nebula as 
a result of interactions of hadrons and leptons (Sect. 4). 
The observed excess of cosmic rays from the direction of the 
Cygnus region at energies of ~ 10^® eV is due to neutrons, 
which are dissolved from heavy nuclei in their interactions 
with the matter of the Cyg OB2 (Sect. 5). We show that 
if a small part of nuclei, which escaped from the PWNa, is 
captured in a strong magnetic field of dense regions of the 
Cyg 0B2 association, the flux of neutrons can explain the 
observed cosmic ray excess from the Cygnus region. 



2 A MODEL FOR HIGH ENERGY 
PROCESSES IN THE CYG OB2 

Let us assume that a very energetic pulsar is formed in a 
core collapse of one of the massive stars in the Cyg 0B2 as- 
sociation. The massive O type stars should evolve through 
the Wolf-Raeyt phase during which the star loses a signifi- 
cant part of its envelope, rotates fast, and has strong surface 
magnetic fields. The Wolf Raeyt stars are probably progen- 
itors of neutron stars with extreme parameters (period, sur- 
face magnetic field) which allow acceleration of particles to 
very high energies. Following the recent works (e.g. Gallant 
& Arons 1994, Blasi, Epstein & Olinto 2000) we assume 
that such neutron stars can accelerate heavy nuclei in their 



wind zones (e.g. in the mechanism called magnetic sling- 
shot, Gunn & Ostriker 1969). Arons and collaborators (see 
e.g. Arons 1998) argue that the iron nuclei can reach ener- 
gies equal to a significant part, Xi of the total potential drop 
through the pulsar polar cap region, 

Bpe = xZe^open, (1) 

where, Ze is the charge of the iron nuclei, ^open = •\/ -Brot /c 
is the maximum potential drop across the polar cup of the 
pulsar, _Erot is the rate of rotational energy lost by the pul- 
sar, X < 1> and c is the velocity of light. These nuclei are 
accelerated with efficiency corresponding to the Goldreich 
& Julian current at the light cylinder radius (Goldreich & 
Julian 1969). Due to the rotational energy losses of the pul- 
sar, the energies of nuclei injected into the expanding nebula 
drop with time. Based on the modelling of the Crab Neb- 
ula, Arons (1998) argues that x should be not very far from 
one. Assuming x = 0.5, we use the following proscription for 
the evolution of the Lorentz factor of injected Fe nuclei with 
time 

7Fe = EFe/rUFe « 4 X 10^ BiiP's , (2) 

where mpe is the mass of the iron rmcloi, B = 10^^-Bi2 G is 
the surface magnetic field of the pulsar, and P = 10~^Pms s 
is the pulsar period which evolves in time t (in seconds) due 
to the rotational energy losses according to 

pL = Po.ms + 2 X 10-^tS?2, (3) 

where Po,ms is an initial period of the pulsar. 

In order to obtain the equilibrium spectra of these nuclei 
inside the nebula at a specific time, we construct a simple 
time dependent model for the expanding nebula, basing on 
the early papers which discuss the evolution of a supernova 
envelope under the influence of a young pulsar (e.g. Ostriker 
& Gunn 1971). This model takes into account possible ac- 
celeration of the nebula due to the energy supply by the 
pulsar. This has influence on the main parameters of the 
nebula, i.e. an expansion velocity, a shock radius, an outer 
radius, and a magnetic field inside the nebula. The details 
of such a model are described in Sect. 2 of Bednarek & Bar- 
tosik (2003). The equilibrium spectra of nuclei arc obtained 
by integrating their injection spectra at a speciflc time over 
the activity period of the pulsar, and taking into account dif- 
ferent processes such as: (1) inelastic collisions of nuclei with 
the matter (important at an early phase); (2) the adiabatic 
energy losses of particles due to expansion of the nebula; 
and (3) the escape of nuclei from the nebula. We include all 
these processes as it has been already done in our previous 
papers (see sect. 4 in Bednarek & Protheroe 2002). However, 
in order not to complicate the model too much, the gravi- 
tational energy losses of the pulsar are neglected. We also 
include partial disintegration of iron nuclei, extracted from 
the neutron star surface, during their propagation through 
the radiation field of the pulsar outer gaps in the inner pul- 
sar magnetosphere following Bednarek & Protheroe (1997). 
Even for the pulsars with the age 2 x 10* years, the iron 
nuclei still lose on average one nucleon in collisions with 
the nonthermal radiation in the pulsar's outer gap. Nucle- 
ons (protons, neutrons) dissolved from the iron rmclci in this 
process are added to the spectrum of nuclei injected into the 
nebula. The example equilibrium spectrum of nuclei, which 
are still inside the nebula at 2 x 10* years after supernova 
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Figure 1. The equilibrium spectra of Icptoiis and different types 
of nuclei (dotted curve) inside the nebula at the time 2 x 10^ yrs 
after supernova explosion. The initial period of the pulsar is 2 ms 
and its surface magnetic field 6 x 10^^ G, and the initial expan- 
sion velocity of the supernova remnant is 2000 km , and its 
mass is 4 Mq . The spectra of Icptons arc shown for two radiation 
fields inside the nebula: (1) the synchrotron photons (SYN) and 
the microwave background radiation (MBR) (full curve), and (2) 
SYN, MBR, and the infrared photons with temperature 100 K 
and energy density 5 times higher than the MBR (dashed curve) . 



explosion, is shown in Fig. 1. A higher energy bump in the 
spectrum of hadrons contains heavy nuclei of the iron group. 
Their energies are hmited to a relatively narrow range due 
to the adiabatic energy losses of nuclei injected at the ear- 
lier phase and the escape of more energetic nuclei from the 
nebula. A lower energy bump in the spectrum of hadrons 
contains protons from disintegrations of heavy nuclei during 
their propagation in the pulsar outer gap and from colli- 
sions of rmclci with a dense matter of the nebula at an early 
stage after supernova explosion. All the medium mass nuclei 
with relatively large Lorentz factors, produced in collisions 
of nuclei with the matter of the nebula at an early phase af- 
ter the explosion, have already managed to escape from the 
nebula. The equilibrium spectra of nuclei inside the nebula 
at a specific time allow us to calculate the spectra of 7-rays 
from decay of pions, produced by nuclei in collisions with 
the matter of the nebula. 

According to Hoshino et al. (1992), relativistic heavy 
nuclei can transfer a part of their energy, ^, to relativistic 
positrons in the process of resonant scattering in the re- 
gion above the pulsar wind shock. As a result, positrons 
are accelerated with the power law spectrum with the in- 
dex a ~ 2 between mc(?^Fc and 'y-pe'm-peC? / Z , where rric 
is the mass of the electron. The injection spectra of these 
positrons evolve in time due to the change of energies and 
injection rate of iron nuclei with the pulsar period. The equi- 
librium spectra of positrons inside the nebula at a specific 
time are calculated by taking into account their radiation 
energy losses on bremsstrahlung, synchrotron, and inverse 
Compton (IC) processes, and due to the adiabatic losses 
caused by an expansion of the nebula as it was shown in 
Sect. 4.2 in Bodnarek & Bartosik (2003). The efficiencies of 
specific energy loss processes are determined by the evolv- 
ing in time parameters of the nebula. The bremsstrahlung 
process depends on the average density of matter inside the 



nebula which can be easily evaluated for the known total 
mass of the nebula and its radius at a specific time. The syn- 
chrotron losses of leptons are determined by the magnetic 
field strength inside the nebula which is estimated from the 
known dependence of the magnetic field inside the pulsar 
inner magnetosphere and the pulsar wind zone (see for de- 
tails Eq. 7 in sect. 2 in Bednarek & Bartosik 2003). The 
evolution of the magnetic field strength in the pulsar wind 
postshock region is shown in Fig. 2b, for the example param- 
eters of the pulsar and the nebula. The IC losses of leptons 
is more difficult to estimate. For relatively old nebulae the 
energy density of background radiation is dominated by the 
microwave background radiation (MBR), and possibly in- 
frared background and optical background. Since the optical 
background created by the massive stars inside the Cyg OB2 
becomes important only for low energy leptons, < 100 GeV, 
we do not consider their influence on the equilibrium spec- 
trum of leptons inside the nebula. The radiation field inside 
relatively young nebulae is dominated by the synchrotron 
photons produced by leptons inside the nebula which spec- 
trum is determined by the conditions inside the nebula in 
the past. Therefore, in order to estimate the density of syn- 
chrotron photons at a specific time we apply simple model 
for the evolution of nebula and include the above defined 
prescription for the injection of particles and the evolution 
of the radiation processes. It is assumed that only a part 
of the whole nebula radiate efficiently. This might be moti- 
vated by the observed anisotropic injection of particles into 
the pulsar wind nebulae which manifests itself by the pres- 
ence of X-ray torus in the equatorial plane of the pulsars 
(e.g. the X-ray torus around ~ 20° of the equatorial plane 
observed in the Crab Nebula). We estimate the density of 
synchrotron photons inside the nebula assuming that lep- 
tons are injected and radiate only within the part of the 
nebula limited by the above mentioned 20° torus. The evo- 
lution of the synchrotron background inside the nebula at a 
specific time after supernova explosion is shown in Fig. 2c. 
It is clear that inside the nebulae older than ~ 3 x 10^ yrs, 
the IC energy losses of leptons on synchrotron photons be- 
comes negligible. Only scattering of the MBR and possibly 
the infrared radiation, produced by the hot gas inside the 
nebula, plays important role. Fig. 2a shows the equilibrium 
spectra of leptons inside the nebula at a specific time af- 
ter supernova explosion. The lower energy part of the spec- 
trum of leptons is formed by leptons which were injected 
relatively early after supernova explosion and lose energy 
at latter time mainly due to the adiabatic expansion of the 
nebula. Therefore, the low energy cut-off in the spectrum of 
leptons, present between 0.1 - 0.01 GeV (Fig. 2a), shifts to 
lower energies with the evolution of the nebula. The efficient 
synchrotron cooling determines the shape of the spectrum 
in the higher energy range when the nebula is young and 
the magnetic field is strong. This process is responsible for 
the power low spectrum with the index close to 2 (seen in 
Fig. 2a) at the time 10^ yrs. On the other hand, the spec- 
trum of leptons at higher energy part at the latter phase of 
expansion of the nebula is formed by freshly injected leptons 
which are not able to cool to low energies fast enough due 
to the weak magnetic and photon fields inside the nebula. 

Described above pulsar-supernova model is applied to a 
variety of high energy processes observed from the direction 
of the Cyg OB2 association. To describe self-consistently all 
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the observed radiation processes, we postulate that the pul- 
sar was born with the initial period Pq = 2 ms and the 
surface magnetic field B = 6 x 10^^ about t = 2 x 10* yrs 
ago, within or close to the dense core of the 0B2 association. 
If this pulsar loses energy mainly on emission of electromag- 
netic radiation, then its present period should be ^ 210 ms 
(see Eq. 3). A pulsar with such parameters resembles the 
Vela type pulsars, i.e. it should produce the pulsed 7-rays 
with luminosity comparable to the Vela pulsar itself or e.g. 
PSR 1706-44. Note that these pulsars are surrounded by 
the synchrotron pulsar wind nebulae from which the steady 
TeV 7-ray emission is observed. It is assumed that the su- 
pernova exploded close to the dense regions of the Cyg OB2 
in the medium with the density of at least ~ 30 cm'^ (Butt 
et al. 2003), and the initial velocity of the bulk mass of the 
supernova envelope is 2000 km s~^ and its mass is 4 M© (as 
derived for the Crab Nebula, Davidson & Fesen 1985). For 
these parameters of the pulsar and the supernova, the ve- 
locity of the supernova front wave after 2 x lO** yrs is ~ 100 
km s~^. Such a low velocity front wave can be difficult to 
observe due to its distortions in the interactions with the 
clumpy high density medium of the Cyg 0B2 association. 
In fact, the supernova remnant shell-like structure has been 
revealed near the TeV source region in the CO, HI, and ra- 
dio emission maps (Butt et al. 2003). The average density 
of matter inside this shell (~ 30 cm~^) is much smaller than 
the average density of matter in the core of the Cyg 0B2 
equal to ^ 300 cm^'', calculated from the total mass of the 
Cyg OB2 ~ 10^ M© (Butt et al. 2003) and its core radius 
~ 15 pc (Knodlseder 2000). 

In the forthcoming sections we argue that the GeV and 
TeV 7-ray emission from the Cygrms region can be due to 
the Vela type pulsar and its nebula, and the excess of ~ 10^* 
eV cosmic rays is due to the neutrons from disintegration of 
very relativistic nuclei which were accelerated soon after the 
pulsar formation and arc still present inside the high density 
region of the Cyg OB2 association. 



3 3EG J2033-h4118 - A VELA TYPE PULSAR 

The spectrum of 3EG J2033-I-4118 was originally described 
as consistent with a simple power law with the index of 
1.96 ± 0.10, but a double power law or a power law with 
the exponential cut-off give better fits (Bertsch et al. 2000, 
Reimer & Bertsch 2001). This source is classified as being a 
non- variable (see discussion in Butt et al. 2003), in contrary 
to its earlier classification as a variable by Mukherjee et 
al. (2003). 3EG J2033+4118 is believed to be located within 
the Cyg OB2 association (Aharonian et al. 2002, Mukherjee 
et al. 2003, Butt et al. 2003). 

Interestingly, the 7-ray emission of this source resem- 
bles the spectra of the middle aged pulsars (the break at 
~ 1 GeV and the sharp cut-off at a few tens of GcV, re- 
quired in the spectrum of 3EG J2033-I-4118 by the level of 
TeV emission observed by the HEGRA telescope). There- 
fore, we suggest that the 7-ray source 3EG J2033-I-4118 is 
caused by a yet undiscovered Vela type pulsar within the 
Cyg 0B2 association. In fact, also the level of the 7-ray 
emission from 3EG J2033-I-4118 is very close to that ob- 
served from the Vela type pulsars. Fig. 3 shows the spectra 
of the Vela pulsar (V), after its re- normalization from the 



Vela pulsar distance ^ 500 pc to the distance to the Cyg 
0B2 ~ 1.7 kpc, and the PSR 1706-44 (P), which is at a 
similar distance as the Cyg OB2. The spectrum of the Vela 
pulsar is measured by the EGRET (Thompson et al. 1997), 
COMPTEL (Schonfelder et al. 2000), OSSE (Strickman et 
al. 1996), and RXTE (Harding et al. 2002). The spectrum 
of PSR 1706-44 is measured by the EGRET (Thompson et 
al. 1996) and constrained by the upper limits on the pulsed 
emission in the X-rays by the ROSAT and ASCA (Becker et 
al. 1995, Finley et al. 1998). It is clear that the 7-ray spec- 
tra of these sources are very similar. The present parameters 
of the pulsar inside the Cyg 0B2 can be like those ones in 
the example considered in the previous section, i.e. the age 
of 2 X 10* yrs, the present pulsar period of ^ 210 ms, and 
the surface magnetic field of 6 x 10^^ G. We encourage the 
searches of such order of periodicity in the 7-ray data of the 
source 3EG J2033-I-4118. 



4 TEV J2032-I-4130 - A PULSAR WIND 
NEBULA 

As it has already been mentioned, the Vela type pulsars are 
surrounded by the non-thermal nebulae observed in the X- 
rays and TeV 7-rays. The present parameters of the pulsar, 
which we propose as respossible for the high energy processes 
in the Cyg OB2, are the closest to the pulsar PSR 1706-44. It 
lies at a similar distance, equal to ~ 1.8 kpc, has the present 
period of 102 ms, and the surface magnetic field of 3.1 x 
10^2 G. The level of non-thermal X-ray emission from the 
nebula around PSR 1706-44, observed by the ROSAT and 
ASCA (Becker, Brazier & Triimper 1995, Finley et al. 1998), 
is close to the level of diffuse X-ray emission observed by 
the Chandra from the direction of the 7-ray source TeV 
J2032-h4130 in the Cyg 0B2 (Butt et al. 2003). 

The X-ray to TeV 7-ray emission from the nebula 
around PSR 1706-44 has recently been interpreted in terms 
of the hadronic, leptonic model for the high energy radiation 
from the pulsar wind nebulae (Bednarek & Bartosik 2003). 
The relatively low level of X-ray emission from this nebula 
(as compared to e.g. the Crab Nebula) is explained if the 
efficiency of positron acceleration by the nuclei injected by 
the pulsar is relatively low ^ w 0.08. On the other hand, 
in order to explain the level of the TeV 7-ray emission, an 
additional target of infrared photons inside the nebula is re- 
quired, apart from the low energy synchrotron photons and 
the MBR (see Fig. 8 in Bednarek & Bartosik 2003). The 
7-rays, produced in collisions of hadrons with the matter of 
the nebula, do not contribute essentially to the TeV 7-ray 
spectrum observed from this nebula due to the relatively low 
density of matter around PSR 1706-44. 

In contrast to the PSR 1706-44 and its nebula, which 
were modeled by the pulsar with the initial period of 15 
ms, here wo postulate that the pulsar in the Cyg OB2 was 
born in a recent supernova explosion of one of the massive 
stars with a much shorter initial period, of the order of 2 
ms. This supernova exploded in a much denser surrounding, 
estimated at ^ 30 cm~^ by Butt et al. (2003), than that 
one considered for the progenitor of the PSR 1706-44. Note 
that the average density of the Cyg OB2 association is sig- 
nificantly higher than the above mentioned value. This can 
be explained by the sweeping effect of the expanding strong 
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Figure 3. The broad band spectrum of the Cyg OB2 association: the diffuse X-ray emission (Chandra - Butt et al. 2003), the GeV 
7-ray emission from the EGRET source 3EG J2033+4118 (Hartman et al. 1999, Lamb & Macomb 1997)), the TeV emission from the 
HEGRA source (Aharonian et al. 2002). The cosmic ray flux at EeV energies is marked by the filled circles (CR). The spectra of Vela 
type pulsars i.e., the Vela pulsar (V) and PSR 1706-44 (P), are shown by the thin dotted curves. The synchrotron and IC spectrum from 
the nebula around the pulsar are shown by the full and dashed curves, respectively. These spectra are reported for two targets of soft 
photons inside the nebula: (1) only the MBR and synchrotron photons (thin dashed curve); (2) additional infrared photons from the 
nebula (thick dashed) . The spectrum of 7-rays from collisions of hadrons with the matter inside the nebula is shown by the dot-dashed 
curve. The spectra of neutrons (included their decay during propagation from the distance of 1.7 kpc) and 7-rays, produced in hadronic 
collisions by nuclei captured in dense regions of the Cyg OB2, are shown by the dot-dot-dot-dashed and thick dotted curves, respectively. 



wind created by the massive, pre-supcrnova star. Because of 
the shorter initial period of the pulsar and the higher den- 
sity medium, the nebulae around the PSR 1706-44 and the 
pulsar in the Cyg 0B2 should differ in some aspects. 

In terms of the hadronic-leptonic model discussed by 
Bednarek & Bartosik (2003), we calculate the synchrotron 
and the IC spectra produced by leptons, which gain energy 
in the resonant interactions with nuclei, and the 7-ray spec- 
tra from decay of pions, produced in the interactions of nu- 
clei with the matter inside the nebula. An example equilib- 
rium spectra of leptons and nuclei inside the nebula 2 x 10'' 
yrs after supernova explosion are shown in Fig. 1. In these 
calculations it is assumed that the surface magnetic field 
of the pulsar born in the Cyg 0B2 association is equal to 
6 X 10^^ G, and its initial period of 2 ms. Note, that the rota- 
tional energy lost by such a pulsar contributes significantly 
to the total energy budget of the expanding nebula. 

The nuclei with the equilibrium spectrum interact with 
the matter of the PWNa producing 7-rays via decay of 
pions. These 7-ray spectra are calculated by applying the 
scaling break model for hadronic interactions proposed by 
Wdowczyk & Wolfendale (1987). The 7-ray spectra emitted 
2 X 10* yrs after supernova explosion are shown in Fig. 3. 
Note that the level of this 7-ray emission is consistent with 
the fiux reported by the HEGRA group. The situation with 
the synchrotron and IC spectra produced by leptons is more 
complicated since the low energy photon target inside the 
nebula is not precisely known. First, we consider the case 
when the soft photon field is created only by the microwave 
background radiation (MBR) (the energy density of syn- 
chrotron photons is already negligible for a 2 x 10* years 



old nebula). The results of calculations of the synchrotron 
and IC spectra from the nebula with the age of 2 x 10* yrs 
are shown in Fig. 3, assuming that the efiiciency for lep- 
ton acceleration by the nuclei is equal to ^ = 0.05. Such 
low efficiency is required in order not to exceed the limit 
on the difl^use X-ray emission from the Cyg OB2 observed 
by the Chandra telescope (Butt et al. 2003). Next, we in- 
crease the soft photon field inside the nebula by adding the 
infrared photons produced by the gas inside or around the 
nebula. The temperature of this radiation is 100 K and the 
average energy density is 5 times higher than that of the 
MBR. The synchrotron and IC spectra calculated with the 
presence of this additional infrared target are shown by the 
thick dashed curve in Fig. 3, using this same efficiency of lep- 
ton acceleration as before. It is evident that only the model 
with the additional infrared target inside the nebula can 
produce IC 7-ray flux on the level consistent with reports 
by the HEGRA group. However, this emission has a much 
steeper spectrum above 1 TeV, spectral index ~ 2.5, than 
that one from the decay of pions. Therefore, the observations 
of TeV J2032-F4130 at energies below 1 TeV should decide 
which processes, hadronic or leptonic, play the main role in 
this source. The comparable contribution of both processes 
to the 7-ray fiux at ~ 1 TeV requires fine tuning of the 
acceleration efficiency of leptons by the nuclei to the value 
mentioned above. 



© RAS, MNRAS 000, 000-000 



6 W. Bednarek 



5 CYGNUS CR EXCESS - A REMNANT OF 
VERY YOUNG PULSAR 

The pulsar with the initial parameters: Po = 2 ms and 
B = 6 X 10^^ G, can generate electric potentials in which 
in principle nuclei can be accelerated above ^ 10^* eV per 
nucloon (see Eqs. 1 and 2). Nuclei with such energies will 
soon leave the nebula around the pulsar and propagate in 
the surrounding medium. The model discussed above for the 
evolution of the PWNa allows us to calculate the radius of 
the nebula and the magnetic field inside the nebula at an 
arbitrary time (see Fig. 2b and also Fig. 1 in Bednarek & 
Bartosik 2003). Nuclei injected inside the nebula at the time, 
iacc, escape from it at the time, tesc, if their diffusion dis- 
tance, Rdif, in the magnetic field of the nebula is equal to 
the dimension of the nebula, -RNcb, at the time iesc- When 
estimating the escape time of rmclei from the nebula we 
take into account the structure of the magnetic field inside 
the nebula (its dependence on the radius), described by the 
model of Kennel & Coroniti (1984). Also the time evolu- 
tion of the magnetic field at the pulsar wind shock is taken 
into account. The diffusion of nuclei through the nebula is 
determined by the diffusion coefficient 

Ddif = (c7^i/3)^/^ (4) 

where Rl = EFe/eZB{t) is the Larmor radius of nuclei, 
Bit) is the magnetic field inside the nebula at the time t. 
The distance travelled by the nuclei is then obtained by 
integration over the time from the moment of injection up 
to the escape (see for details Bednarek & Protheore 2002), 

^dif = \l^dt'. (5) 

The above mentioned condition, ildif = ^^Neb, allows us to 
determine the moment when nuclei with different energies 
and injection time escape from the nebula. The spectra of 
nuclei which managed to escape from the nebula during the 
first 2 X 10"* yrs after explosion of supernova are shown in 
Fig. 4, assuming the initial parameters of the pulsar and 
nebula as considered above. Our calculations show that nu- 
clei with the Lorentz factors above --^ 10*' should diffuse out 
of the nebula into the surrounding medium during the first 
2 X 10* yrs. 

We consider the situation in which the pulsar (and its 
nebula) is formed inside or close to the high density regions 
of the Cyg 0B2 association. The nuclei with the highest en- 
ergies, Lorentz factors ~ 10^, which escaped into the galactic 
medium with the magnetic field of ~ 5 x 10~^ G, have the 
Larmor radii of ~ 200 pc, and diffuse during 2 x 10* yrs up 
to R^if ~ 600 pc from the Cyg 0B2. It is likely that a part 
of these nuclei can be captured by dense regions for ~ 10* 
years. The magnetic fields in dense molecular clouds are 
much stronger than in the Galactic medium, and can reach 
the values of ~ 10"* - 10"^ G (e.g. Crutcher 1999). In such 
magnetic fields, the Larmor radius and the diffusion distance 
of nuclei with the Lorentz factors ~ lO'* are ~ 10 pc and 
~ 140 pc during the time of 2 x 10* years. These distances 
are close to the radius of the central core in the Cyg 0B2 
estimated on Rob2 ^ 15 pc and to its diameter of the equal 
to ~ 60 pc. Therefore, since the pulsar is close to the Cyg 
OB2, we estimate that up to ^ ~ {RoB2/R^i^)~^ « 10"^ of 
all injected nuclei can be captured by the nearby Cyg OB2. 



The spectra of nuclei which escaped from the PWNa 
are shown in Fig. 4. Nuclei with different mass numbers 
were produced in fragmentation of heavy nuclei (the iron 
group), occurring during their collisions with the matter of 
the expanding supernova remnant during the first several 
years after explosion, when the column density of matter 
was still large. Note that two groups of heavy nuclei can be 
distingtiished in the total spectrum (see Fig. 4). The higher 
energy heavy nuclei escape from the supernova envelope just 
after their injection into the nebula by the pulsar since their 
Lorentz factors are large enough that corresponding Lar- 
mor radii are comparable to the dimension of the nebula at 
that time. On the other hand, the lower energy heavy nu- 
clei escape from the nebula at a later time (~ 10'* years), 
after suffering significant adiabatic energy losses due to the 
expansion of the nebula. 

We calculate the flux of neutrons produced in collisions 
of nuclei with the matter inside the Cyg OB2 association. 
The average density of matter in the Cyg OB2 is estimated 
at ~ 300 cm""*, for the core radius of the Cyg 0B2 equal 
to 15 pc and its total mass of ~ 10® Mq. It is assumed 
that in a single collision, the nuclei are split into two lower 
mass stable nuclei with the rest of nucleons released free. 
The energies of these nucleons are reduced by a factor of 
two in comparison to energies of nucleons in primary nuclei. 
The flux of neutrons produced in this process is shown in 
Fig. 3. Interestingly, the neutron flux calculated in terms of 
the model discussed here, is consistent with the excess of 
cosmic ray events at energies ~ 10^* eV, reported from the 
Cygnus region by the AGASA group equal to ~ 8% of the 
cosmic ray flux (Hayashida et al. 1999b). 



6 CONCLUSION 

We have shown that different types of high energy phenom- 
ena observed from the direction of the Cyg 0B2 association 
can be self-consistently explained by the model of energetic 
pulsar created in the recent explosion of one of the massive 
stars. The pulsar should have the initial period close to 2 ms. 
For the present age of 2 x 10* yrs, its period should be close 
to ~ 210 ms, if its surface magnetic field is 6 x 10^^ G. Such 
a pulsar belongs to the family of the so-called Vela type pul- 
sars, producing pulsed 7-rays with the intensity comparable 
to the EGRET source, 3EG J2033-H4118, observed in the 
direction of the Cyg 0B2. In fact, the spectral features of 
3EG J2033-I-4118 (possible break in the spectrum, required 
high energy cut-off) are characteristic for the spectra of the 
Vela pulsar and the PSR 1706-44. Therefore, we encourage 
the searches of the hundred millisecond periodicity in the 
7-ray sources reported in the direction of the Cyg OB2. 

The pulsar in the Cyg 0B2 should also be surrounded 
by the pulsar wind nebula, as observed in the ceise of other 
Vela type pulsars. We propose that the TeV 7-ray source re- 
ported by the HEGRA group is due to the 7-rays produced 
by leptons and/or hadrons accelerated by the PWNa in- 
side or close to the Cyg OB2. Therefore, in spite of the close 
physical relation, the TeV 7-ray emission from the Cyg OB2 
(produced in the PWNa) do not need to connect smoothly 
with the GeV 7-ray emission (produced in the pulsar mag- 
netosphere) . 

We also suggest that the cosmic ray excess, observed 
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from the Cygnus region by the AGASA experiment, is also 
physically related to this energetic pulsar. The highest en- 
ergy nuclei, injected by the pulsar soon after its formation, 
escape from the PWNa. A part of these nuclei can be cap- 
tured by the strong magnetic fields in the high density re- 
gions of the Cyg OB2. These nuclei collide with the matter of 
the Cyg OB2, injecting neutrons which create the observed 
excess of cosmic rays in the direction of the Cyg OB2. 

The 7-rays from the Cyg 0B2, produced in collisions 
of hadrons with the matter, should also be accompanied by 
the comparable fluxes of neutrinos. However, the predicted 
muon neutrino event rate in a 1 km^ detector of the Ice- 
Cubc type during a 1 yr observation is rather low (based on 
the calculations of the likelihood of detecting neutrinos by 
such a detector obtained by Gaisser & Grillo (1987), and the 
absorption coefficients of neutrinos in the Earth derived by 
Gandhi (2000)). We estimate that hadrons captured in the 
dense regions of the Cyg OB2 should give between ~ 0.14 
to ~ 0.05 neutrino events per year, for the source observed 
close to the horizon (not absorbed by the Earth) and from 
the Nadir direction (absorbed by the Eaxth), respectively. 
The hadrons confined in the PWNa should produce ~ 0.48 
(~ 0.46) events per yr for the source located at the horizon 
(the Nadir). The spectra of neutrinos produced by hadrons 
in the PWNa and the Cyg OB2 are shown in Fig. 5. Neutrino 
spectra produced in the Cyg 0B2 are above the atmospheric 
neutrino background expected within 1° of the source but 
they are below the sensitivity of the 1 km'^ neutrino detector. 
On the other hand, the neutrino spectra produced inside the 
PWNa are closer to the sensitivity limit of the 1 km^ detec- 
tor but they are significantly below the atmospheric neutrino 
background. Therefore, we do not predict any detection of 
neutrinos from the direction of the Cyg 0B2 association. 

Finally, we conclude that the presence of at least one 
Vela type pulsar in the Cyg 0B2 association is quite likely 
on statistical grounds, since the supernova rate in this asso- 
ciation should be one per a few tens of hundreds years on 
average (the Cyg OB2 contains about 100 O type stars with 
the lifetime of a few milion years) . 
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Figure 2. (a) The equilibrium spectra of leptons inside tlie neb- 
ula, injected with the power law spectrum as discussed in the 
main text, calculated for the pulsar with an initial period of 2 ms 
and a surface magnetic field of 6 X 10^^ G, and the supernova 
with a mass of 4Mq, expanding with an initial velocity of 2000 
km s~^, are show for different time after a supernova explosion 
(and a pulsar formation) lO'^ yrs (dot-dashed curve), 3 X 10^ yrs 
(full), 10^ yrs (dashed), 3 X lO^ yrs (dot-dot-dot-dashed), and lO" 
yrs (dotted), (b) The magnetic field strength in the pulsar wind 
postshock region (inside the nebula) as a function of time, (c) 
The differential density of synchrotron photons inside the nebula 
at a specific time (the curves are marked as in figure (a)). 




log (E^uci / GeV) 

Figure 4. The spectrum of rmclei (thick curve) which escaped 
from the pulsar wind nebula into the surrounding medium. 
The spectra of nuclei with difi'crciit mass immbers are shown 
by: the thin full curve (protons), dashed curve (A=2-10), dot- 
dashed curve (A=ll-20), dotted curve (A=21-40), and dot-dot- 
dot-dashed curve (A=41-55). It is assumed that the pulsar was 
born 2 X 10* yrs ago with the initial period of 2 ms and the sur- 
face magnetic field of 6 x 10^^ G, and the initial velocity of the 
supernova remnant is 2000 km s~^ and its mass is 4 Mq. 
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Figure 5. The spectrum of muon neutrinos, produced by hadrons 

captured in the high density regions of the Cyg Ob2 associa- 
tion (full curve) and produced by hadrons inside the PWNa (dot- 
dashed curve), are compared with the atmospheric neutrino back- 
ground (dashed curves) within 1° of the source for the zenith and 
azimuth directions (Lipari 1993). The 3-yr sensitivity of the Ice- 
Cube detector is marked by the dotted line (Hill 2001). 
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